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We study the excitations in a three dimensional version of Kitaev's spin-^ model on the honeycomb 
lattice introduced by the present authors recently. The gapped phase of the system is analyzed using 
a low energy effective Hamiltonian H^ff which is defined on the diamond lattice and consists of 
plaquette operators. The excitations of H^ff form loops in an embedded lattice. The elementary 
excitations, which are the shortest loops, are fermions. Moreover, the excitations obey nontrivial 
braiding rules: when a fermion winds through a loop, the wave function acquires a phase tt. 
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I. INTRODUCTION 

Kitaev's spin-^ model on the honeycomb lattice has 
become a paradigmatic system in the study of topolog- 
ical order in condensed matter systems^. The model is 
exactly solvable; there is a gapped phase that supports 
Abelian anyons and a gapless phase in which the exci- 
tations (in the presence of a perturbing and gap induc- 
ing external magnetic field) are non- Abelian anyons. A 
remarkable feature of Kitaev's Hamiltonian is that, in 
contrast to other topologically ordered solvable models 
such as the toric-code^, it involves only local two-spin 
interactions, which makes it a physically feasible model. 
Various aspects of the model have been extensively stud- 
ied so far—"—. 

The key element in Kitaev's construction, leading to 
its exact solvability, is that the lattice is trivalent and 
that the links can be labeled with three different colors in 
such a way that at each site no two links have the same 
color. Using this fact, the present authors have gener- 
alized his construction to three dimensions'^. There ex- 
ists other three dimensional generalizations of the Kitaev 
modeUii^. 

The 3D Kitaev model introduced in Ref. [13 has also 
been solved exactly and has a phase diagram similar to 
the 2D case, with a gapped phase and a gapless one. 
It has further been shown that the gapped phase is de- 
scribed by an effective toric-code-like Hamiltonian -ffe// 
defined on the diamond lattice and consisting of plaque- 
tte operators. In this paper, we analyze the excitations 
of Heff- We find that they have the structure of loops, 
the energy of a loop being proportional to its length. 

The elementary excitations, which are along the small- 
est loops, are fermions. Since the excitations arc one di- 
mensional objects, it is meaningful to ask, even though 
they live in three dimensions, whether they obey nontriv- 
ial braiding statistics. We find that the state acquires a 
phase TT when a fermion winds through a loop. 

The paper is organized as follows. Section [H] briefly 
reviews the three dimensional Kitaev model and the low 
energy effective Hamiltonian Hgf f describing its gapped 
phase. In Sec. IIIII we give the representation of the exci- 



tations of Heff as loops in an embedded lattice. Calcu- 
lation of exchange and braiding statistics is presented in 
Sec. IIVI We conclude with a discussion of our results in 
Sec. El 



II. HAMILTONIAN 




FIG. 1: The 3D lattice: The four sites inside the loop (marked 
1 to 4) constitute a unit cell; ai , a.2 and as are the basis 
vectors. Plaquette p consists of sites marked 1 to 10. 



The Hamiltonian is defined on the lattice shown in Fig. 
[1] For details about the lattice we refer to Ref. Here 
we just note that the coordination number is three and 
the links, which connect neighboring sites, are labeled 
X, y or z. In addition, the three links at each site have 
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different labels. The Hamiltonian is 

<ij>x <i-j>y <i,j>z 

(1) 

where trf are the Pauli matrices at site i and J2<ij> 
denotes summation over a-typc links. 

The above Hamiltonian can be solved exactlj*i^ using 
a Majorana fermion representation of Pauli matrices^. 
The excitation spectrum is gapless when Jz Jx + Jy, 
Jx < Jy + Jz, and Jy < Jz + Jx- Everywhere else in the 
parameter space there is a gap. The purpose of this paper 
is to study the nature of low-energy excitations in the 
gapped phase. Specifically, we will concentrate on the 
limit Jz > Jx,Jy 

Let -ffo be the Hamiltonian obtained by putting 
Jx ~ Jy ^ in Eq. ([T]); it corresponds to isolated z- 
links and the spectrum is trivially solved. The ground 
state has a large degeneracy: any state in which the two 
spins in any given z-link, say, connecting sites i and j, 
are either | tt)ij or I has the minimum energy. Here 
t)8 = I t)», <| l)j = - | and similarly for j. 

For small values of Jx and Jy, the low energy excita- 
tions (more precisely, those eigenstates of H adiabatically 
evolving from the degenerate ground states of Hq, as Jx 
and Jy are turned on) are described by an effective Hamil- 
tonian Hef f acting on this degenerate subspace (denoted 
by Ho). Such an effective Hamiltonian will be defined on 
the lattice formed by the z-links, which we call Ceff, 
and it turns out to be the diamond lattice, as far as the 
connectivity is concerned (see Fig. [2]). Ceff, which has 
coordination number four, is obtained from the original 
lattice by shrinking each z-link to its mid-point. 

The links in Cef / lie along four different directions. We 
can choose an orthogonal coordinate system such that 
these directions are (1,0,1), (—1,0,1), (0,-1,-1) and 
(0,1,-1). We label the links respectively as a,b,c and 
d. (Note that the links are not directed; positive and 
negative directions are equivalent.) Then at each site, 
the four links have four different labels. 

Let / denote the sites of L^ff- Ho is spanned by the 
states { H/ l"^/)7 | m/ e {+, -}}, where |+)/ = | tt)/ 
and = I 4-4-)/. The local degrees of freedom acting 
on Ho are effective spin-i operators rf (a £ {x,y,z}) 
defined by 

t;|+)/ = |+)/, r;|-), = -|-)z, 

rf\+)i = \-)i, r^\-)i = -\+)i- (2) 

Before writing down H^ff, it is convenient to define 
a loop operator as follows. Consider a loop I in Cef / of 
length n and consisting of sites {Ik\k = 1, 2, ... n} such 
that Ik and Ik+i are nearest neighbors, with the identifi- 
cation In+i = Ii- Let 7fe denote the link connecting sites 
Ik and Ik+i- We then define 

Bi = f[rT:, (3) 

fe=i 




FIG. 2: The lattice Ceff (diamond lattice) on wliicli the 
effective Hamiltonian is defined. Bold, long-dashed, medium- 
dashed and short-dashed lines respectively denote links with 
labels a, b, c and d. 



where ak is determined by the labels of the links "fk-i 
and jk as follows: 

ak — X for labels (a, c) and (6, d), 

ak — y for {a,d) and (6, c), 

ak = z for (a, 6) and (c, d). (4) 

Here we need to consider only pairs with different labels 
since "fk-i and jk are both links involving site k and all 
links at a given site have different labels. It is easy to see 
that, with this definition, for any two loops / and /', 

[Bi,Br]=0. (5) 

This can be understood as follows. First, note from Eqs. 
(|4]) that ak is the same for two pairs of labels if and 
only if they are complementary to each other. Now, I 
and I' can intersect each other at one or more sites, and 
they can also have overlapping segments. Let / be a site 
where they intersect. Then the two pairs of labels for 
the links that meet at /. corresponding to I and I' , will 
be complementary to each other and therefore ak will be 
the same for both the loops. Now consider a segment 
where they overlap. Again, ak will be the same for I 
and I' everywhere on the segment except at the two end 
points, where they are strictly different. This is true for 
any overlapping segment of I and I' , which in turn means 
that among those sites shared by the two loops, there 
will be an even number of them where ak is different for 
the two loops. This implies Eq. ([S]) since rfrf — —TjTf 
when a ^ b. 

To lowest order in Jx and Jy , Hef / is given in terms of 
the above operators corresponding to the smallest loops 
in Ceff - We call such loops plaquettes and use the index 
p to denote them. A plaquette consists of six sites and 
there are four types, each type corresponding to a par- 
ticular label sequence and orientation. The four possible 
sequences arc acbacb (type-1), cadcad (type-2), bdabda 
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(type-3) and dbcdbc (type-4). In Fig. [2 {1,2,3,4,5,6}, 
{7, 8, 9, 3, 2, 1}, {8, 10, 5, 4, 3, 9} and {7, 1, 6, 5, 10, 8} are 
respective examples. 

In terms of the plaquette operators, the lowest order 
effective Hamiltonian isi^ 

^ " 25^7^ W'-^y ^ + J"4T.'bp) : (6) 

^ \ p p / 

where the unprimed sum is over plaquettes of type- 1 and 
2, while the primed sum is over plaquettes of type- 3 and 
4. 



III. EXCITATIONS 

Since all plaquette operators commute among them- 
selves [see Eq. ([5|)], i?e// is, in principle, trivially diago- 
nalized. The energy eigenstates are then labeled by the 
eigenvalues of Bp, which are -1-1 or —1. However, not all 
Bp arc independent. In fact, it is the constraints existing 
among them that make Heff nontrivial and give rise to 
topological excitations. 

Next we write down all the constraints and then find 
the allowed configurations of Bp that are consistent with 
the constraints. In the process we will also calculate the 
ground state degeneracy. From now on we assume peri- 
odic boundary conditions in all three directions. 

We first find a complete set of mutually commuting 
conserved operators that includes a maximal number of 
plaquette operators. To this end, we define a prod- 
uct of two loops. For our purposes, it is convenient to 
represent a loop by the set of links it contains. Let 
k = {71,72, •■ -Tni} and I2 = {<5i, (^2, • . • ^nj be two 
loops in Ceff containing ni links (denoted 7^) and 712 
links (denoted 6i) respectively. The product of li and I2 
is then defined as 

li2 = li-l2 = liUl2-linh. (7) 

That is, li2 is the loop formed by the links contained in 
h or I2 but without those links common to both. Then 
it is easy to check that 

Bh2 = Bi^ ■ Bi^. (8) 

Therefore, to form a complete set of conserved operators 
we need to consider only the elementary plaquettes, but 
not any loop that can be obtained by combining a sub- 
set of the latter. There are three topologically nontrivial 
(noncontractiblc) loops corresponding to the three direc- 
tions, which cannot be obtained from the plaquettes; we 
denote them Cx , C'y and Cz respectively. Any loop in the 
lattice can be obtained by combining the plaquettes and 
the three noncontractiblc loops. However, as we stated 
earlier, the plaquettes themselves are not all indepen- 
dent. There are three types of constraints that plaquette 
operators satisfy. 



a. Local constraints. It follows from Eq. ([8|) that 
the product of any set of plaquettes that forms a closed 
surface equals 1. The smallest such surface is formed 
by four adjacent plaquettes, any two of which share two 
links, as shown in Fig. [2] We call the volume enclosed by 
such a surface a basic cell. There are TV basic cells, where 

is the number of sites. This gives rise to — 1 inde- 
pendent constraints; it is A^ — 1 because the constraint 
corresponding to any one basic cell equals the product of 
the constraints corresponding to the remaining A^ — 1. 

b. Surface constraints. As a consequence of Eq. ([8]) 
there are three further constraints, independent of the 
local constraints, which correspond to the three topolog- 
ically nontrivial surfaces. 

c. Volume constraints. Any product of corner shar- 
ing plaquettes that fill out the whole lattice equals one. 
There is only one such independent constraint. 

Putting everything together, the total number of con- 
straints is A^ -I- 3. There are 2A^ plaquettes and the num- 
ber of independent plaquette operators is A^ — 3. These 
along with the three noncontractiblc loop operators Cx, 
Cy and Cz form a complete set of A commuting opera- 
tors. 

There is another set of conserved operators correspond- 
ing to noncontractiblc surfaces. Such an operator can be 
defined on a plane S that divides the lattice without cut- 
ting any links in such a way that at each site lying on 
the plane, two of the links lie on one side of the plane 
and the remaining two on the other. Then the surface 
operator is defined as follows: 

Ms = n^r, (9) 

les 

where the product is over sites lying on S and a/ is deter- 
mined by the two incoming (or outgoing) links according 
to the rule given in Eq. Q. Mg commutes with all con- 
tractible loop operators because noncommuting — in our 
case, anticommuting — contributions come only at points 
where the loop pierces the surface and there will be an 
even number of such points. There are three such op- 
erators corresponding to the planes normal to x, y and 
z direction. We denote them Mx, My and M^, respec- 
tively. 

Just as we obtained loop operators as a product of link- 
sharing plaquette operators, we can construct conserved 
'surface' operators by multiplying corner sharing plaque- 
tte operators. In such a product, contribution from a 
shared site cancel out since it will be the same Pauli ma- 
trix for the two plaquettes sharing it [owing to Eq. 
As a result the product will only involve sites that are not 
shared. Consequently, when corner sharing plaquettes fill 
out a certain volume in the lattice, the contribution to 
the operator will come only from the sites on the bound- 
ary. 

The set of noncontractiblc loop and surface operators 
Cx: Cy, Cz and Mx, My, Mz form a system of three 
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spin-i degrees of freedom since 

CaMa = -MaCa, CaAh = M^Ca for a ^ b, 

CaCb = CbCa, MaMb = A'hMa. (10) 

It then follows that the ground state degeneracy is 2^ — 
8; since the Hamiltonian does not depend on the above 
nonlocal operators. In fact, all excitations will have the 
same degeneracy of 8. 

In the gromid state \GS), Bp ~ +1 for all plaquettes. 
Evidently, all the constraints are satisfied. The excita- 
tions are obtained by flipping some of the Bp's to — 1. 
To find the configurations of Bp that are consistent with 
the above constraints, it is useful to consider the lattice 
formed by the basic cells, each of which correspond to 
four plaquettes forming a closed surface. This lattice, 
denoted Cp, is also the diamond lattice and the links in 
Cp represent the plaquettes in the original lattice. 

In Cp, the local constraints have a simple picture: the 
plaquette operators corresponding to the four links at 
each vertex multiply to 1. This implies that at every ver- 
tex only an even number of links (0, 2 or 4) can be fiipped 
to -1, which in turn means that links with Bp = — 1 form 
closed loops in Cp. It is easy to see that any such loop 
excitation will also satisfy the surface constraints. How- 
ever, the volume constraint imposes a further restriction 
on the allowed loop excitations. 

The smallest loop in Cp is an elementary plaquette 
consisting of six links; we call it a 6-loop to distinguish 
it from the plaquette in Ceff- Nevertheless, the exci- 
tation of an odd number of 6-loop violates the volume 
constraint. This can be shown by explicit construction, 
but there is another more transparent way to see this. 

We have seen that the loop operators in Eq. ^ leave 
the ground state invariant. However, an open-string op- 
erator, obtained from a loop operator by truncation, i.e. 
by removing the part corresponding to a segment, is not 
conserved. The noncommuting part is at the ends of the 
string; it anticommutes with 6 of the Bp's involving the 
spin at a given end. The plaquettes corresponding to 
these 6 Bp's form a 6-loop in Cp. An open-string oper- 
ator then creates two 6-loop excitations at the ends. If 
two excited loops li and I2 overlap, then the resultant ex- 
citation will be along the loop I12, which is the product 
of li and I2 [see Eq. ([7])]. This is because the overlapping 
links are flipped twice — equivalent to no flip. Since any 
noncontractible loop / can be expressed as a product of 
6-loops, the excitation along / can be created by exciting 
the corresponding 6-loops. 

In the extreme case in which the string consists only of 
one site, say /, the corresponding string operator is rf , 
rf or rf . Then the two 6-loops at the 'ends' overlap (two 
of the links are shared), and the resultant excitation is 
an 8-loop. Since any state in the Hilbert space must be 
obtained from the ground state by the action of a combi- 
nation of r" operators, it follows that all states will have 
an even number of 6-loops, the elementary excitations. 

However, a single 6-loop can be isolated by taking the 
other end of the string that creates it far away. This does 



not cost energy since the loops do not interact except 
when they overlap (then the overlapping part is annihi- 
lated). The energy of a loop is proportional to its length, 
which corresponds to the number of Bp's excited. 

In other words, the excitation spectrum has two super- 
selection sectors: 1) the vacuum, which consists of states 
with an even number of 6-loop excitations and 2) the 
6-loop, which consists of states with an odd number of 
6-loop excitations. 

In general, an excitation corresponds to a configuration 
of one or more loops in Cp, where the loops can intersect 
but not overlap. As we have just shown, such an exci- 
tation is created by exciting the 6-loops into which the 
loop configuration can be decomposed. Let i be a loop 
configuration and let F be one of the decompositions of L 
into 6-loops. It is straightforward to explicitly construct 
the operator that creates the excitation along L: arbi- 
trarily pair up the constituent 6-loops; connect each pair 
by a string, again arbitrarily; then take the product of 
all the string operators that excite each pair. 

Here it should be stressed that only the loop configu- 
ration L, on which the excitation lives, is physical. The 
state depends neither on the particular choice of decom- 
position r nor on the specific way in which the 6-loops 
are paired and connected by strings. All choices give rise 
to the same state, up to a phase. 

A loop excitation can also be interpreted as the edge of 
a membrane. Just as open-string operators are obtained 
by truncating a loop operator, we can construct an open- 
membrane operator by removing a patch from a closed- 
surface operator defined as a product of corner sharing 
plaquettes. Let Om be such an operator defined on a 
membrane M. As in the case of string operators, Om 
does not commute with the Hamiltonian. Nonetheless, 
it is only certain Bp's along the edge of M that do not 
commute with Om] in fact, they anticommute. There- 
fore, in the state Om\GS) these anticommuting Bp's are 
excited; they will correspond to a loop in £p. In the lan- 
guage of the previous construction, Om corresponds to a 
specific choice of decomposition into 6-loops, which fixes 
M, and then neighboring 6-loops are paired together and 
connected by the shortest string between them, which 
consists of a single site. 

IV. STATISTICS 

In this section, we determine the statistics obeyed by 
the excitations. We will see that the elementary 6-loop 
excitations are fermionic. Nontrivial phases also arise 
from braiding: when a 6-loop winds through a bigger 
loop, the state acquires a phase tt. 

A. Localized excitations 

A 6-loop excitation is the maximally localized state 
permitted by the constraints. So, they can be treated 
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as particle-like and then one can ask whether they are 
fermions or bosons. We recall that in three dimensions 
there is no other statistics for point particles. 

In a lattice, the statistics is determined by the algebra 
of the hopping operatorsi^. Let i,j,k,l be four sites in 
the lattice and tais, a, (3 £ J, fc, be the correspond- 
ing hopping operators. Then, 

iikiiih] = iijiuhk for bosons, 

iikhiii] = -iijiidik for fermions. (11) 

A 6-loop in Cp represents 6 plaquettes in Ceff, the 
lattice on which H^ff is defined. In Ceff, such a set 
of 6 plaquettes share a common link. Moreover, there is 
no other plaquette of which this particular link is part of. 
This provides a unique one-to-one mapping between links 
in Ceff and 6-loops in Cp. Therefore 6-loop excitations 
can be thought of as localized around the links in Cef / . 

The above mapping gives a description of loops in Cp 
in terms of a link variable in Cef / : it takes value — 1 if 
the corresponding 6-loop is excited, otherwise, -1-1. This 
description, however, is redundant due to a Z2 gauge 
symmetry related to the non-uniqueness of the decompo- 
sition of a loop into 6-loops. The gauge transformation 
at a given site flips all the link variables connected to 
that site. 




FIG. 3: A site and the four links connected to it in the lattice 
Ceff 



Consider a site / in Cef / and the four links emanating 
from it, denoted i,j,k,l (see Fig. [3]). In the previous 
section, we showed that rf excites an 8-loop, which can 
be decomposed into two 6-loops. By the above mapping, 
this corresponds to a pair of links at /, the value of a 
determining which pair. Note that, due to the gauge 
symmetry, exciting a given pair or its complement are 
equivalent. 

If r" excites the links, say, i and j, then it also hops 
an excitation between the two links. Without any loss of 
generality, we can assume that the links are labeled such 
that the hopping operators are 



tij 


= tkl 




tik 


= iji 


= rl 


til 


= ijk 





The first equality in each of the above equations follows 



from the gauge symmetry. Moreover, iy- ~ iji. It imme- 
diately follows that 

tiktiitij = —tijtiitik- (13) 

Therefore, 6-loops arc fermions. 

B. Braiding 

A configuration of loops can evolve along topologically 
distinct paths. The question one can ask then is whether 
there is a phase difference between distinct paths. We 
will now show that a phase tt is acquired when a 6-loop 
winds through a bigger loop. 

Let Om be a membrane operator that creates an exci- 
tation along loop I. 

leM 

Here M is a set of sites / which fills the surface enclosing 
I, and a I is determined by the labels of the two incoming 
(or outgoing) links, according to Eq. (|4]). 




FIG. 4: Schematic depiction of a fermionic 6-loop braiding 
through a bigger loop. Traversing such a path results in a 
topological phase tt 

A 6-loop is moved along a closed path by the loop 
operator corresponding to the path. A path which takes 
it through / will have one site common with M. Here 
we take / to be big enough so that the 6-loop can pass 
through without interacting, i.e., without touching (see 
Fig. 13]). Let this site be Ic- The loop operator Bi is given 
by 

lei 

where /?/ are also determined by the set of rules in Eq. 
©• 

Since / goes through the membrane, /S/^ is determined 
by one outgoing and one incoming link; therefore, it is 
strictly different from a/^ . Hence 

OmBi = -BiOm- (16) 

This results in a phase difference of tt between paths 
which go through I and paths which do not. 
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V. SUMMARY AND DISCUSSION 

We have studied the gapped phase of 3D Kitaev model 
by analyzing the effective Hamiltonian H^f f obtained in 
the limit of strong z— links. Hf.ff is defined on the di- 
amond lattice and the individual terms it contains arc 
mutually commuting plaquette operators. There are loop 
and surface operators that are conserved. In the language 
of Levin and Weni^, our model has both string-net and 
membrane-net condensation. 

There is an embedded lattice Cp, also a diamond lat- 
tice, in which the links stand for plaquettes. In Cp, the 
excitations have the structure of loops. The energy of 
a loop is proportional to its length, consequently, the 
ground state is free of loops. The loops interact only 



when they overlap, in which case the overlapping part 
is annihilated. The smallest loops in Cp, the 6- loops, 
are the elementary excitations and they obey fermionic 
statistics. The fermionic 6- loop braids nontrivially with 
a bigger loop: when it winds through the latter over a 
closed path, the wave function acquires a phase tt. 

Even though there already exists other exactly solvable 
models in 3D that are topologically ordere d^'^i^^" — , our 
model has a particular relevance in that it is obtained 
as the effective Hamiltonian for the gapped phase of 3D 
Kitaev model, which involves only two-spin interactions 
as opposed to higher-spin operators in the other solvable 
models. It will therefore be worthwhile to analyze the 
model at finite temperature to see if topological order 
survivesii— . 
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